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Abstract 

During cosmological inflation, it has been suggested that fields coupled 
to the inflaton can be excited by the slow-rolling inflaton into a quasi- 
stable non-vacuum state. Within this scenario of "warm inflation", this 
could allow for a smooth transition to a radiation dominated Universe 
without a separate reheating stage and a modification of the slow roll 
evolution, as the heat-bath backreacts on the inflaton through friction. In 
order to study this from first principles, we investigate the dynamics of a 
scalar field coupled to the inflaton and N light scalar boson fields, using 
the 2PI-l/iV expansion for nonequilibrium quantum fields. As a first step 
we restrict ourselves to Minkowski spacetime, interpret the inflaton as a 
time-dependent background, and use vacuum initial conditions. We find 
that the dominant effect is particle creation at late stages of the evolution 
due to the effective time-dependent mass. The further transfer of energy to 
the light degrees of freedom and subsequent equilibration only occurs after 
the end of inflation. As a consequence, the adiabatic constraint, which is 
assumed in most studies of warm inflation, is not satisfied when starting 
from an initial vacuum state. 



^email: g.aarts@swan.ac.uk 

§ email : A . Tr anberg@damtp . cam. ac . uk 

1 



1 Introduction 



In recent years quantitative cosmological observations have become available, re- 
sulting in a model for our Universe that is consistent with an early period of 
inflation and predictions of simple inflationary theories [I]. In the most elemen- 
tary setup, the dynamics of the inflaton mean field <p(t) = (<p(t, x)) is determined 
by 

i(t)+3H(t)i(t) + V'[<f>(t)] = 0, (1.1) 

subject to the Friedmann equation H 2 = (^<fi 2 + V r [0])/3Mp 1 , where Mpi is the 
Planck mass. 

In warm inflation [21 El SI El El [T| the key idea is that interactions between 
the inflaton and other quantum fields are important during inflation and that 
they result in continuous energy transfer from the inflaton to these other fields. 
If this transfer is sufficiently fast and equilibration is rapid, a quasi-stable state 
could be achieved, different from the inflationary vacuum. An additional effective 
friction term would then be expected in the inflaton equation of motion. A simple 
phenomenological modification of Eq. (11.11) incorporating this reads [2], [3] 

4>(t) + [3H{t) + T>(t)] 0(t) + V'[<j>(t)} = 0, (1.2) 

where T^(t) is a time- and field- dependent friction coefficient FPI 

Since warm inflation necessarily involves multiple interacting quantum fields, 
which are dynamically evolving in real time, a full, quantitative understanding is 
difficult. In particular, a first-principle investigation requires all tools available 
to study quantum field dynamics far from equilibrium beyond the mean-field 
approximation. This may be contrasted with the theory of preheating due to 
parametric resonance after inflation, which can be understood from a combination 
of mean- field dynamics and the classical random field approximation [TU | [Til [T2"] . 

For the purpose of this paper, the scenario of warm inflation separates natu- 
rally into three parts [5]: 

1. The inflaton field tp is coupled to a second scalar field x- While the inflaton 
rolls down the effective potential, it interacts with the x field, resulting in 
the excitation of x degrees of freedom. 

2. The field x is coupled to other degrees of freedom, which can be light 
fermion (?/>) or scalar (a) fields. These fields are excited as well and may 
thermalize. Various interaction terms are possible. 

3. The backreaction of the x ne ld on the inflaton leads to effective friction in 
the inflaton equation of motion, resulting in overdamped inflaton dynamics. 

1 For a critical analysis of warm inflation, sec Rcf. [8]. 

2 For a study of particle creation and friction during power law inflation, see Ref. [9j. 
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A commonly used interaction term between the inflaton field and the second 
field x is given by ^g 2 tp 2 x 2 - Due to the time dependence of the inflaton (f)(t), 
an effective time-dependent mass term appears in the quantum dynamics of the 
X field, mt + g 2 <j) 2 (t), where m x is the mass in absence of the inflaton. In this 
paper, we consider a large field model, such that the inflaton is rolling down its 
effective potential and <f>(t) is decreasing. As a first step, we treat the inflaton as 
a time- dependent background and assume the most extreme case, in which the 
slow-rolling of the inflaton is caused by the interactions with the heat-bath rather 
than the expansion of the Universe. Taking this overdamped limit with T 3> H, 
and assuming a quadratic potential V[ip] = \m} ip Lp 1 for simplicity, the dominant 
solution of Eq. (jl.2p in the overdamped regime reads <p{t) = <f>o exp[— (m 2 /T)t\. 
We assumed that T is time- and field-independent. The inflaton is treated as 
a dynamical field in a companion paper [13]. Here, our goal is to study the 
dynamics of parts 1.) and 2.) described above, by combining a mode function 
analysis [2] and the techniques of the two-particle irreducible (2PI) effective 
action [15] for nonequilibrium quantum fields [16]. Since the inflaton acts as a 
time-dependent background, we find that there is particle production in the x 
sector, akin to particle production in curved spacetime [H]. In Section [2] we 
study this process using a mode function analysis. We find that details of x 
particle creation depend crucially on the size of the zero-temperature mass m x 
and the momentum k. Most importantly, we find that particle production only 
takes place towards the end of the evolution and that the amount of particles is, 
to a large extent, independent from the initial conditions and the duration of the 
inflationary stage. 

In order to determine the range of validity of the mode function analysis 
and study the effect of interactions, we continue in Section [3] with a full far- 
from-equilibrium numerical study in quantum field theory using the 2PI effective 
action, and include interactions between the scalar field x an d N light quan- 
tum fields. Specifically, we include N scalar fields a a (a = 1, . . . , N) and use a 
truncation of the 2PI effective action determined by the 2PI-l/iV expansion to 
next-to-leading order (NLO) [T71 ITS]. This approximation has been well stud- 
ied in recent years and is known to give a quantitative description of both the 
early evolution far from equilibrium as well as the later stages of equilibration 
and thermalization (see e.g. Refs. [TTJ, [191 EH [21]). We consider both a trilinear 
coupling x a a an d a quartic coupling x 2(T l- Irrespective of the interaction term, 
we find that as long as the inflaton is evolving, the mode function analysis gives 
an accurate description of the dynamics in the x sector. In particular, processes 
leading to equilibration and thermalization are not yet relevant. 

In this first study, we ignore the expansion of the Universe and work in 
Minkowski spacetime. If anything, neglecting the dilution caused by expansion 
should increase the chances of realising a "warm" state. 
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2 Mode function analysis 

As a first step, we perform a mode function analysis for a free \ field, subject to 
a time-dependent mass 

M 2 (t) = m\ + 5mV 7 *. (2.1) 

Here 5m contains the details on the coupling to the inflaton. Specifically for the 
\g 2 tp 2 "X 2 interaction, we find that 5m 2 e -7 * = g 2 <f) 2 (t), but these details are not 
required here. Using the standard decomposition, we write 

d 3 k 



X(t, x) 



ak/k(t)e ikx + 4/*(t)e 



t t*i 



kx 



(2tt) 3 

and find that fk(t) is a solution of 

f k (t) + [k 2 + m 2 + 5m 2 e" 7t ] f k (t) = 0. 



(2.2) 



(2.3) 



A change of variables to x = (25m/ 7 */ 2 shows that this is a Bessel equation 
of the form 



2 AuJ V 

x H £ 



/_ 



/k(x) = 0, 



(2.4) 



where u; k = ^/k 2 + m 2 . We note that only the combination ur^/'y appears in this 
equation, while the dependence on 5m will enter via the initial conditions. The 
general solution is 



25m 
e 

7 



-7*/2 



25m 
< 

7 



-7*/2 



(2.5) 



where J^z) is the Bessel function of the first kind. The constants are de- 
termined by the initial conditions, which are fixed by demanding that /k(0) = 

l/\/2T\ and / k (0) = -ifi k /k(0), where fi k = ^k 2 + M 2 (0) and M x (0) is the 

initial mass. This determines as 



A ± 
A k 



where 



TV 



y/2Tt^ 7 sinh (27ru; k / 7) 



J: 



T2io; k /7 



25m 

7 



k^q=2io; k /7 



25m 

7 



J: 



=F2icj k /7 



25m 

7 



dt J ^h 



25m 
< 

7 



-7*/2 



, (2.6) 



(2.7) 



t=o 



One may verify that the Wronskian f]i(t) f^(t) — fk(t) f^_(t) = i is preserved during 
the time evolution. 

In order to identify the produced particle number at asymptotically late times, 
we use some elementary properties of Bessel functions and find that 



lim J. 

t— >oo 



±2iw k /7 



25m 
( 

7 



- 7 t/2 



°k e ' 



(2-« 




Figure 1: Particle number as a function of according to the Bogoliubov 
relation for four values of the final mass m x . The initial mass is M x (0)/7 = 20. 



with 



5m\ 
7 / 



1 



r(l ± 2za; k / 7 ) (2 ' 9) 

At late times, we find therefore that the mode functions oscillate with the ex- 
pected frequency u^, and 



lim / k (t) = A+C+e'^ + A-C-e^. 



t— >oo 



(2.10) 



A comparison with the standard form of the mode functions at t — > 00, /k(£) 



/ V2u;k yields the Bogoliubov coefficient [H 



Mt)d t Mt) - Mt)dJ k (t) = -V2^a-c } 



k ■ 



(2.11) 



Therefore, the final particle number (starting in vacuum initially) is n k = |/?k| 2 
2a>k I ArC^ 1 2 , which, after some algebra, can be written aa^| 



x k w k 
1 



7T 



2fik 7 sinh(27TCi;k/ 7) 



J- 



2iirijJk /■y 



25m 

7 



J: 



25m 

7 



1 

~ 2" 
(2.12) 

The resulting particle numbers are shown in Fig. [1] for four different values of 
the asymptotic mass m x and an initial mass M x (0)/7 = 20. We observe that a 
significant amount of particles are only produced when m x / r y <C 1, and only with 
momentum fc/7 < 1. 



3 If the initial particle number is nonzero and equals (a k <2k) = n£\ we find that 

;i + 2n k °V k |c k -| 2 (l^l 2 + l^l 2 )-i 
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Figure 2: Particle number as a function of time, according to the Bogoliubov 
relation, for three values of the initial mass M x (0) and three momentum modes. 
After shifting time by an amount 2 log[M x (0)/M x ef ], where M£ ef /7 = 20 is a 
reference mass, there is no dependence on the initial mass remaining. The final 
mass is rn x /^ = 0.1. 



In order to determine when particles are created during the inflationary stage, 
we show the time-dependent particle number, 

n k (t) = ^7-hA(t)| 2 + ^(t)|/ k (t)| 2 l -I, (2.13) 



2ufc(t) 
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where u^it) = yk 2 + M 2 (t), in Fig.[2j We comment on possible other definitions 

of particle number elsewhere [13]. The time-dependent particle number is shown 
for three different initial masses M x (0), the largest mass corresponding to the 
longest period of inflation (recall that M x (0) ~ </> )- However, we find that a 
trivial shift of the time variable is sufficient to take the initial mass dependence 
into account. In other words, the amount of particles produced is independent of 
the initial conditions and therefore the duration of the inflationary stage. Gener- 
ically, particles can only be produced when pk(^)| ^ a; k(^); which in our model 
translates to 7 > M x (t). If we denote the time when most particles are created 
with we find from Fig. that - 2 log[M x (0)/M x cf ] » 8, where M x ef / 7 = 20 
is a reference mass. The size of the effective mass is then M x (t*)/^y ~ 0.37 < 1, 
confirming that \ particles are only produced when M x (t)/~f is sufficiently small. 



3 Interactions 

In the mode function analysis of the previous section, interactions besides the 
time-dependent mass are ignored. In order to assess the validity of that analysis, 
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we now include interactions and consider the coupling of the scalar field x to 
TV light degrees of freedom, either fermionic or bosonic. We start with some 
parametric estimates and then present the numerical results obtained with the 
help of the 2PI-l/iV expansion for nonequilibrium quantum fields. 



3.1 Parametric estimates 

In order to get an estimate for the magnitude of various parameters, we use 
here an inflaton potential V[ip] = ^m^ip 2 . During slow- roll inflation, the Hubble 
parameter is approximately given by H ~ m^/Mpi. In order to have inflation, 
the initial amplitude should be <po ~ a fewMpi, and to satisfy CMB constraints 
P, we require that m v ~ KT 6 M P1 . We then find that H ~ lO" 6 0o > 10 -6 M P i. 
In warm inflation, the effective damping term T in the inflaton equation wins over 
the expansion rate. Taking T ~ 100H yields T ~ 10 _4 M P i. Using exponential 
time dependence, </>(t) = 0o exp[— (m 2 /T)t\, then yields a value of 7 ~ 10~ 8 M P i 
for the rate in the time-dependent mass (12. ip . In order to have any particle 
production, we found from the mode function analysis that m x /j < 0.1, or 
m x < 10 _9 Mpi, which means that the x particle can be some degree of freedom 
beyond the Standard Model. 

We now consider the couplings to the light degrees of freedom. First we 
consider the coupling to N fermion fields, with the interaction term 

N , 

^2-7fiX&a4>a- (3.1) 
a=l * 

The factor \fN is introduced to allow for a proper 1/N expansion. Following 
Ref. [5], we impose h 2 < 1 such that perturbation theory is reasonable!^ From 
a standard one-loop calculation one can compute the (zero-temperature) onshell 
decay width for x ~ > m the case that M x > 2M^. Here the masses M Xj ^ 
include possible background field dependence. The width is given by 

h 2 M 2 f 4M|\ 3/2 
^ " Sir^TJq V " ~Ml) ■ (3 ' 2) 

The important assumption made in most studies of warm inflation is that the 
dynamics takes place in the so-called adiabatic approximation, \4>N>\ T£, lead- 
ing to quick decay (and possibly thermalization) during inflation^ In our model 
we have to compare T£ with 7. Taking for simplicity <C M x and p = 0, we 



find 



YJ = I^M X 

7 8n 7 



(3.3) 



4 This corresponds to Nh 2 < 1 in the conventions of Ref. [S]. 

5 The origin of the adiabatic approximation can be traced back to Refs. [22l [23] 
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When it is assumed that M x ~ g<po and 0o 3> 7 (as discussed above), this ratio 
can be much larger than one. However, in the mode function analysis we found 
that x particles are only produced when the (time-dependent) mass M x is much 
smaller, specifically M x /j < 1. In that case the important conclusion is that the 
dynamics is not taking place in the adiabatic regime, but rather in the opposite 

limit 10/01 > r* 

In the numerical study below, we couple the x field to N light bosonic fields, 
with the interaction term 

i=l * 

In this case, the coupling constant h is dimensionful and it is natural to write 
h = m x h, where h is dimensionless. We consider again the decay process x ~ > a<J 
and find for the onshell width, assuming that M x > 2M a , 



h 2 m 2 x / 4M 2 

vl = ^Wrw x Y~^t' (3 ' 5) 

To test the adiabatic approximation, we compare again T£ with 7. Taking M a <C 
M x and p = 0, we find 

T% h 2 m x m x 

— = « ~M~ « L ( 3 - 6 ) 

7 Sir 7 M x 

Since all factors are strictly less than one, we are not in the adiabatic limit. This 
is in agreement with the results from the numerical analysis carried out below. 



3.2 Nonequilibrium dynamics 

To determine the range of validity of the mode function analysis and estimates 
carried out above, we continue with a numerical study using the 2PI effective ac- 
tion for quantum field dynamics in real time. We refrain from providing details on 
the 2PI effective action and the Schwinger-Keldysh formalism for nonequilibrium 
field theory, instead we refer to Ref. |18j . whose notation we follow closely. 
We consider the following action, 



s \x> a \ = ~ 




We consider two different interaction terms Mnt[X;°"] between the x fi e ld and 
the scalar fields a a , a trilinear and a quartic coupling. Both interaction terms 
preserve the O(N) symmetry in the a sector, such that the a two-point function 
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can be written as G aab (x,y) = 5 a bG a (x,y). We have scaled all couplings in such 
a way that a proper 1/N expansion is possible. We use N = 4 throughout. 

Trilinear coupling. Motivated by Ref. [5], we start with the trilinear coupling, 

Vtoblx, <r] = ~7^X^l + c xX- (3-8) 



Because this term breaks the symmetry X ~Xi we have to allow for a nonzero 
expectation value x(t) = vN (x(t, x))Jj The term linear in x is used to shift the 
minimum of the (effective) potential at the initial time to (x) = 0. The 2PI part 
of the effective action is written as T 2 [G X , G a } = T xa [G x , G a ] + r a [G a ], where the 
first term is given by the two-loop diagram 

T xa [G x ,G a ]=ih 2 j d 4 xd 4 yG x (x,y)G 2 a (x,y). (3.9) 

The second part, r 2 [G CT ], is the standard NLO contribution for N scalar fields. 
This contribution has been discussed in detail in Refs. [T71 [TS] and will not be 
shown here explicitly. Other 2PI diagrams are suppressed by 1/N and are not 
included. The resulting equations of motion follow from a variation of the 2PI 
effective action with respect to x, G x and G a , and read 



d 2 

dT^ + 



X = — hG a (x, x) — c x , (3.10) 



and for the two-point functions (j = x, °~) 

- p x + Mf(t)]G j (x,y) = i J d 4 zZ,(x,z)G J (z,y) + i5 4 (x-y), (3.11) 
with the effective masses 



M 2 x (t) = m x + We" 11 + -fx, M&t) = m 2 a + 2h X + X^-^-G^x, x). 

(3.12) 

The nonlocal contributions to the self energies S X . CT follow from variation of 
T 2 [G X , G a ] in the usual manner [TH [T3] . 

Quartic coupling. In order to assess the importance of the possibility of on- 
shell decay x ~^ ao ~ m th e trilinear case, we also consider the following quartic 
potential, 

v*AxA = ^°l (3-13) 

In the 2PI effective action we include the two- and three-loop diagrams 

h f ih 2 f 

T xa [G x ,G a ] = -- / d 4 xG x (x,x)G a (x,x) + — / d 4 xd 4 y G?(x, y)G 2 a (x, y), 

(3.14) 



^This was overlooked in Ref. [5]. 
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Figure 3: Equal-time two-point function G x (t, t; k) as a function of time t in 
units of 7, for nine different momenta k (the ordering is such that the smallest 
momentum corresponds to the largest value at ft ~ 10), using the analytical ex- 
pression (top left), the trilinear interaction (bottom left), the quartic interaction 
term in the Hartree approximation (top right), and the quartic interaction term 
(bottom right). In all cases the initial mass is M x (0)/7 = 20 and the final mass 
is rn x j^i = 0.1. 



as well as IV [G^] to NLO as above. Strictly speaking, we deviate here from 
the 1/N expansion, since the three-loop diagram only appears at next-to-next-to 
leading order. We include it nevertheless, since it plays the same role as the 
two-loop diagram in the trilinear case and results in \o interaction beyond the 
mean-field approximation. In this case the symmetry \ ~^ ~X is preserved, so 
that we can take (%) = consistently for the entire evolution. Hence Eq. (13.101) 
is absent, while the effective masses appearing in Eq. (13.111) are 

M*(t) = m\ + 5mV 7t + hG a (x, x), (3.15) 

Mlit) =m 2 a + \o-—^j^—G a {x, x) + ^G x (x, x). (3.16) 

The nonlocal self energies follow again from the 2PI effective action. 

We solve the resulting equations numerically, following the approach in [TTJ 
[T9t [20], [2Tj . Space is discretized on a three-dimensional lattice with 32 3 sites and a 
physical size of 7L = 32. The length of the memory kernel is ft = 20, containing 
800 time steps. The masses are m x / r y = 0.1 and m a /m x = 1/3. The coupling 
constants are X x = 6, A CT = 6. In the trilinear case the coupling h/m x = 5/3, 
while in the quartic case we use h = 1. We initialize the two-point functions in 
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Figure 4: Same as in Fig. [3j for the a two-point function G a (t, t; k). The largest 
value at t — corresponds to the smallest momentum, since G a (0, 0; k) = |(k 2 + 
ml)' 1 . 



vacuum. Renormalization is carried out in such a way that the set of equations 
is initialized at the fixed point of the renormalized mean field equations [2U ITS] . 

Particle number is a derived concept and not always well defined in an in- 
teracting field theory out of equilibrium. Instead of comparing time-dependent 
particle numbers, we prefer to study basic quantities appearing in the dynamical 
equations: the equal-time two-point functions G x (t, i;k) and G a (t,t; k). In Figs. 
[3] and H] we show the evolution of the nine lowest momentum modes in time (the 
evolution of the zero momentum mode is not shown). In the top left corners, the 
dynamics without interactions is presented, in which case G x (t,t;k) = |/k(£)| 2 
and G a (t,t; k) = |(k 2 + m 2 ,) -1 . The latter is exactly time-independent, because 
of our choice of initialization. In the other three frames, we show the dynam- 
ics with trilinear interaction (bottom left), and with quartic interaction using 
the Hartree approximation, i.e. including the self-consistently determined masses 
(I3.15P and (13. 16[) but not the nonlocal terms (top right), and the full evolution 
(bottom right). 

As can be seen in Fig. |3j the evolution of G x is in all cases nearly identical to 
the free case, except for a small additional growth in the case that the nonlocal 
self energies are included. From this we conclude that the mode function analysis 
describes the important aspects of this part of the dynamics extremely well, even 
in the presence of interactions. In particular, the process of particle production 
can be studied in a satisfactory manner using the techniques of Section [21 

The response of the a propagator is shown in Fig. HI Recall that without 
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interactions, the equal-time propagator is constant in time. In the trilinear case 
interactions with the x field have a small effect, in the low- momentum modes 
only and with a characteristic time scale much longer than I/7. In particular, 
there is no notion of equilibration and thermalization. This is consistent with the 
estimates given above and also with previous studies of nonequilibrium scalar field 
dynamics using the 2PI effective action [TTI 1201 |2T] , in which it was shown that 
thermalization occurs only after a time of the order of at least several hundred 
elementary oscillations. In the quartic case, the situation is slightly different 
because of the presence of the x tadpole in the effective a mass ( 13. 16f) . During 
the evolution, the mass M a grows as G x (x, x) increases due the time-dependent 
inflaton background. The response of the o propagator is captured very well in 
the Hartree approximation. Additional nonlocal diagrams (Fig. HI bottom right) 
play only a minor role at this stage, since in the theory with a quartic coupling 
onshell decay x ~^ X ' ' is kinematically not allowed, while scattering processes 
XX aa are suppressed due to small amount of x particles present. 

4 Outlook 

In order to investigate the dynamics of warm inflation, we studied the quantum 
evolution of a scalar field x coupled to the inflaton as well as N light scalar 
fields. In this first analysis we treated the inflaton as a background field and 
ignored the expansion of the Universe. We took vacuum initial conditions. From 
a comparison with nonequilibrium quantum field dynamics using the 2PI-l/iV 
expansion to next-to-leading order, we found that the response in the x sector 
can be accurately understood from a mode function analysis and that further 
interactions with light scalar fields are subdominant. The important conclusion 
is that x particle production only takes place towards the end of the evolution, 
independent of the duration of the inflationary stage. An immediate consequence 
is that under these conditions the so-called adiabatic approximation, \<fi/ <fi\ <C T x , 
assumed in most studies of warm inflation, is violated. Therefore, the dynamics 
during inflation takes place far from equilibrium and processes important for 
thermalization become important only later. These findings are complementary 
to those obtained in Ref. [8]. 

Although our investigation used a specific <p 2 x 2 interaction and assumed an 
exponential time dependence for the inflaton background field, we believe that 
most of the results obtained here are generic for systems initially in vacuum. In 
particular, x particles can only be produced when the rate of time variation of 
the effective mass is comparable with the effective mass itself. This observation is 
a potential hurdle in all models where the effective x mass receives a contribution 
from interactions with the inflaton and the inflaton initially has a large expec- 
tation value (large field models). Concerning the time dependence, we have also 
studied a linear evolution, <p{t) ~ t, and found that the specific time dependence 
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is not crucial [T3] . 

As a next step, we plan to treat the inflaton as a dynamical quantum field, 
and use 2PI effective action techniques to analyse the effective inflaton equation 
of motion, including the backreaction from x particles. Since we found in this 
paper that the role of the additional light degrees of freedom is subdominant, we 
can focus entirely on the inflaton- x sector. 
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